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DISSOCIATION  EQUILIBRIUM 

ABSTRACT 

Results  of  calculations  carried  out  for  a  model  of  air  in  dissociation 
equilibrium  are  presented  in  graphical  form.  The  quantities  computed  are 
i)  flow  variables  (including  species  concentrations)  behind  normal  and 
oblique  shock  waves,  ii)  flow  variables  in  axisymmetric  conical  flow 
fields,  iii)  stagnation  point  values  of  flow  variables  on  the  'stagnation' 
streamline  behind  two-dimensional  and  axisymmetric  detached  shock  waves, 
and  iv)  flow  variable  gradients  at  the  shock  wave  on  stagnation  streamlines. 
Computations  are  given  for  free  stream  temperatures  of  273.l6°K  and  300°K, 
free  stream  pressures  of  1.0,  .1,  .01,  .001,  and  .0001  atmospheres,  and  a 
range  of  initial  Mach  numbers  and  cone  angles  to  provide  flow  field 
temperatures  in  the  range  3000°K  -  10,000°K.  Brief  derivations  of  the 
equations  employed  are  given. 

The  present  calculations  are  oriented  toward  application  in  experi¬ 
ments  in  hypersonic  flow  with  ground  facilities  such  as  shock  tubes  and 
ballistic  ranges.  In  addition,  they  furnish  important  supplementary  in¬ 
formation  to  theoretical  studies  of  nonequilibrium  flows. 
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SYMBOLS 


Kjf. 


h 

h-i 


V 


kb. 

J 


K. 

J 

K 

w 

n 


q 

R 

R 

w 


s 

T 

U 

V 

Wi 

W 

00 

B 

y 


Functions  occuring  in  Eqs  (2.7)  and  (2.9) 

Concentration  of  i ^  species  [k  mol/kg  of  mixture] 

Specific  heat  of  air  at  constant  pressure 

[Dyn  m/kg  deg  K]  (l  Dyn  =  1  Newton  =  10^  dynes) 

=  dh^/dt  (see  Table  II )  [Dyn  m/k  mol  deg  K] 

Specific  heat  of  air  at  constant  volume 

[Dyn  m/kg  deg  K] 

Specific  enthalpy  [Dyn  m/kg]  =  ^  ci  hq 

Molar  specific  enthalpy  of  i^*1  species  (see  Table  II ) 

Reaction  rate  constants  for  reaction,  forward  and  backward, 

respectively  (see  Table  i) 

(sk^  /  k^  ).  Equilibrium  constant  for  reaction  (see  Table  i) 

j  0 

Shock  wave  curvature 
Free  stream  Mach  number 

Arc  length  along  curve  normal  to  streamlines  [m,  mm] 

Pressure  j^Dyn/m2J 

Flow  velocity  [m/sec] 

Universal  gas  constant  =  8312.4  [Dyn  m/k  mol  deg  K] 

Shock  wave  radius  of  curvature  (  =  l/l^) 

Radial  polar  coordinate  ( -  [*2  +  A1'2) 

Arc  length  along  streamline  [m,  mm] 

Temperature  [deg  K] 

Velocity  component  in  r  direction  (  =  q  cos  [0  -  cd]) 

Velocity  component  in  cp  direction  (  =  q  sin  [0  -  cp]) 

Molecular  weight  of  i^1  species  [kg/k  mol] 

Molecular  weight  of  air  (  =  28.8587  gm/mol) 

Angle  between  shock  wave  and  x-axis  [radians,  deg] 

Ratio  of  specific  heats  of  air  (cp/c^) 
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e 

e 

p 

CT 

CP 

b 

i 

j 

t 

w 


=  0  for  two-dimensional  flow;  =  1  for  axisymmetric  flow 

Angle  between  streamline  and  x-axis  [radians,  deg] 

Density  ^kg/m^J 

Arc  length  along  shock  wave 

Angular  polar  coordinate  (cp  =  arctan  y/x) 

Subscripts 

Body  surface 
i^  species 
reaction 

Stagnation  point 
Shock  wave 
Free  stream 
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1.  INTRODUCTION 


Current  experimental  research  in  hypersonic  flow  over  two-dimensional 
and  axisymmetric  bodies  makes  it  desirable  to  have  available  calculations 
for  air  in  chemical  (dissociation)  equilibrium.  Primarily  these  are  values 
of  the  flow  variables  (including  concentrations  of  chemical  species)  behind 
shock  waves  and  in  conical  flow  fields.  Furthermore,  it  is  desired  to  have 
these  data  for  ranges  of  free  stream  conditions  applicable  to  experiments 
with  ground  facilities,  such  as  shock  tubes  and  ballistic  ranges.  To  this 
end  extensive  calculations  of  shock  wave  quantities  and  conical  flows  have 
been  carried  out  on  the  BRL  high-speed  computers  for  a  model  of  air  (see 
Section  2)  in  chemical  equilibrium,  and  results  are  presented  graphically 
here  in  Figs.  1.1  through  II. 9. *  Brief  derivations  of  the  equations  are 
given  in  Sections  3  and  4.  Although  computations  of  the  above  have  already 
been  carried  out  by  other  authors  (e.g.,  Refs.  1  through  4  and  the  refer¬ 
ences  contained  in  them),  the  present  paper  furnishes  information  hitherto 
not  available,  to  the  authors'  knowledge,  in  a  form  convenient  for  work  in 
hypersonic  ground  facilities. 

The  present  calculations  also  supply  important  supplementary  informa¬ 
tion  for  the  determination  of  nonequilibrium  dissociating  airflow  over  wedges 

s**  r 

and  cones.  Certain  features  of  these  nonequilibrium  flows  (e.g.,  entropy 

layers,  oblique  shock  equilibrium  regions)  require  knowledge  of  equilibrium 

value  s . 

Studies  being  conducted  at  BRL  on  the  subsonic  region  in  front  of  a 
supersonic  blunt  body  by  analysis  of  interferometric  data  suggest  that  theo¬ 
retically  determined  information  on  the  stagnation  streamline  would  be  use¬ 
ful.  Therefore,  extensive  calculations  were  performed  to  obtain  stagnation 
values  of  the  flow  variables,  and,  in  addition,  gradients  of  the  flow 
variables  along  the  stagnation  streamline  at  the  shock  wave.  These  results 
are  presented  graphically  in  Figs.  III.l  through  IV. 8;  derivations  are  given 
in  Sections  5  and  6. 


*The  computer  program  is  available  for  cases  not  explicitly  graphed  in 
this  report  should  exact  values  be  required. 

1 

**  Superscript  numbers  denote  references  which  may  be  found  on  page  56. 
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2.  AIR  MODEL 


The  air  is  considered  to  be  a  mixture  of  neutral  species  consisting  of 
0  and  N  atoms,  and  Np»  Og,  and  NO  molecules.  This  mixture  of  particles  is 
assumed  to  be  in  translational,  rotational,  and  vibrational  equilibrium  at 
all  times.  Effects  of  electronic  excitation  and  vibration-dissociation 
coupling  are  neglected.  The  chemical  reactions  are  listed  in  Table  I. 

Let  c^  denote  the  concentration  of  the  i^h  species  (moles  of  per 
unit  mass  of  air)  and  K-<  the  equilibrium  constant  of  the  j^*1  reaction.  K. 
can  be  determined  quite  accurately  from  quantum  statistical  calculations, 
and  the  results  are  often  fitted  to  an  equation  of  the  form 


kd  =  AjT 


n  . 
D 


exp(-Ej/T)  , 


where  T  is  temperature,  and  A.,  n.,,  and  E.  are  constants.  The  values  listed 

J  J  J 

for  equilibrium  constants  in  Table  I  are  based  on  data  given  in  Refs.  6  and  7. 


The  law  of  mass  action  for  equilibrium  flow  leads  to  the  relations 


CN0 


(2.1) 


(where  o  is  the  density  of  the  air)  plus  the  following  restrictions: 


L1(T)  = 


.  *6*7 
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K?K7 

=  1,  l2(t)  =  =  1, 


l3(t) 


.  _  , 

-  %  -  1  ' 


Calculation  shows  that  L^>  Lg,  and  L^  =  1  ±  .10  for  a  limited  temperature 
range,  the  results  being  valid  only  within  this  range  (~  4000°  —  8000°K). 

The  flows  studied  here  are  produced  by  objects  moving  at  constant 
supersonic  speed  in  stationary  air,  taken  to  be  a  mixture  of  two  ideal 
gases  0g  and  Ng  having  concentrations 


. 21153  mole 

°2  Wm  S®  air 

00 


_  .78847  mole 

N„  W  gm  air 

00 

00 
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where  the  subscript  «  denotes  free  stream  conditions,  and  W^,  is  the  molecular 
weight  of  air  (=  28.85870). 

The  conservation  of  chemical  elements  leads  to  the  following  relations: 


Substitution  of  Eq. 


(2.l)  into  Eq. 


+  (1/2)  cQ  +  (1/2)  c^ 

+  (1/2)  cN  +  (1/2)  c^  • 

(2.2)  gives 

pc  c 

+  (1/2)  c0  +(1/2) 

5 

+  (1/2)  cN  +  (1/2)  PCNC0  . 

Kr- 


(2.2) 


(2.3) 


Considering  each  component  as  an  ideal  gas,  the  equation  of  state  for 
the  ith  species  is  (p,  p,  and  T  being  pressure,  density,  and  temperature, 
re spectively) 


and  for  the  mixture  (where  p  =  ^  Pi> 


=  RP1T/Wi 


p  =  RpT^  Ci  (with  c^  =  (p./p)/W.)  .  (2.4) 

i 

"fch 

R  is  the  universal  gas  constant,  Wj_  the  molecular  weight  of  the  i  species, 
c^  the  concentration.  By  Eq.  (2.2) 
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(2.5) 


)_  ci  -  (VO  +  (l/2)(c0  ♦  Cj,)  • 
i 

The  differentiated  versions  of  Eqs.  (2.4)  and  (2.3)  give  the  following 
useful  set  of  linear  equations  for  dp,  dc^,  dc^,  and  dT: 


Audp  +  A12dc0  +  A13dcN  +  Al4dT  =  (2/p)  dp  (2.6a) 

A21dp  +  A22dc0  +  A23dcN  +  A24dT  =  0  (2.6b) 

A3ldp  +  A32dco  +  A33dcN  +  A34dT  =  ®  *  (2.6c) 


Expressions  for  the  coefficients  A  are  presented  in  the  Appendix. 

A  fourth  relation,  valid  along  streamline q ,  if  obtained  from  the 

energy  equation 


dh  =  (l/p  )dp  (2.7) 

where  h  is  enthalpy  per  unit  mass.  The  total  enthalpy  is  the  sum  of  the 
enthalpies  of  the  components: 

h  =Zcihi  • 

i 

Expressions  for  h^  (enthalpy  per  mole)  and  dh^/dT 
Table  II  for  all  the  species.  Eq.  (2.7)  becomes 

A4ldP  +  A42dc0  +  A43dcN  +  A44dT  =  (V°)  dP  »  (2.9) 

the  coefficients  appearing  in  the  Appendix. 

Eqs.  (2.l)  are  differentiated  to  give 

dc02  =  (C0/Ki)[c0^  +  2pdc0  -  (pc0/K1)(dK1/dT)dT] 


(2.b) 


(=  c  )  are  given  in 


pi 
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dcN  *  (cl/K2)[CNdD  +  2pdcN  -  (P%/K2)(dK2/dT)dT] 


(2.10) 


dCH0  '  ^  [C0CN  ^  +  p(cNdC0  +  C0dCN)  '  .... 


!Vn  ^  dT 

K,.  dT 


3.  SHOCK  WAVE  CALCULATIONS 


The  conditions  immediately  behind  a  shock  (denoted  by  the  subscript  w) 
are  given  by  the  following  relations  obtained  from  the  conservation  laws* 
(referring  to  Fig.  3.l): 


w 


w 


tan  (B  -  0w)  =  tan  3 

(3.1) 

cos  (3  -  9  )  =  q  cos  3 

W  co 

(3.2) 

=  P»  +  2  (!  -  Pco/P*)311128 

(3.3) 

=  \  +  d/2)^  [1  -  (poo/Pw)2]sin2B 

(3.4) 

where  fl  is  the  angle  of  inclination  of  the  shock  wave,  q  the  flow  speed, 
and  6  the  angle  of  inclination  of  the  flow.  In  the  free  stream 

•>.  =  (7/2)htA  ,  =  BejJv.  (3.5) 

p„%,2/p.  =  'fjfj  . 

where  M  is  the  Mach  number  and  y  is  the  ratio  of  specific  heats,  having 

CO  00 

the  value  1.4. 

The  flow  variables  behind  the  shock  wave  are  determined  by  solving 
Eqs.  (3.1),  . ..,  (3.4),  together  with  the  equation  of  state  Eq.  (2.4),  plus 
Eqs.  (2.l)  and  (2.3).  These  form  a  set  of  ten  functional  equations  for  the 
ten  variables  pw>  V  Tw>  Pw>  <<=o2)w,  (^J,.  (=0>„.  <CN>W. 

when  the  parameters  M,T,p  ,  (orp)  and  3  are  given.  The  system  of 

00  CO  00  CO 


*The  basic  jump  conditions  for  a  steady  oblique  shock  wave  requiring  con¬ 
servation  of  mass,  momentum,  and  energy  can  be  found  in  many  places;  e.g., 
p.  8  of  Ref.  8. 
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equations  is  solved  on  the  BRL  high  speed  computers  by  successive  application 

of  the  method  of  "regula  falsi,"  or  "false  position."  Frequently  0  is 

w 

given,  and  8  is  the  unknown  quantity;  an  additional  iterative  procedure  (e.g., 

regula  falsi)  can  then  find  the  8  corresponding  to  a  given  9  . 

w 

Figs.  1.1  through  I. 11  contain  curves  of  flow  variables  behind  normal 
and  oblique  shock  waves.  Pressure,  temperature ,  density,  and  species 
concentrations  are  conveniently  expressible  as  functions  of  the  parameter 
sin  8  for  given  free  stream  temperature  and  pressure.  For  the  flow  de¬ 
flection  9,  another  parameter,  M  ,  is  required. 

00 

4.  CONICAL  FLOW 

In  axisymmetric  conical  flow  a  straight  cone  of  half  angle  9^  gives  rise 
to  a  straight  attached  shock  wave  inclined  at  angle  8.  It  is  convenient  here 
to  introduce  polar  coordinates  r,  cp  (as  in  Fig.  3.1.b)  and  to  employ  U  and  V, 
the  components  of  velocity  in  the  r  and  cp  directions,  respectively.  The 
values  of  cp  at  the  shock  and  body  are,  respectively 

With  the  condition  that  the  flow  variables  be  independent  of  radius 
(d/dr  =  0),  the  mass  and  momentum  conservation  relations  reduce  to 

dU/dcp  =  V  (4. 2. a) 

dV/dcp  =  -  (v/p)dp/dcp  -  [2U  +  V  cot  cp]  (4.2.b) 

dp/dcp  =  V2  dp/dcp  +  pV[U  +  V  cot  cp]  (4.2.c) 

On  substituting  dp/dcp  from  Eq.  (4.2.c)  into  the  right  hand  sides  of 
Eqs.  (2.6)  and  (2.9)  one  obtains  four  linear  algebraic  equations  for  dp/dcp, 
d-Co/^Pj  dc^/dcp,  611(1  dT/dcp,  which  are  solved  to  give  differential  equations 

dp/dcp  =  F1?  dcQ/dcp  =  F2,  dcN/dcp  =  F^  dT/dcp  =  F^  (4.3) 

where  F^  F2,  F^,  and  F^  are  functions  of  cp,  U,  V,  p,  p,  T,  cQ,  and  c^. 

Eqs.  (4.2),  (4.3),  and  (2,10  )  form  a  set  of  ten  first  order  differential 
equations  for  U,  V,  p,  p,  T,  cQ^,  c^,  cQ,  cN,  and  c^,  which  can  be 

l4 
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integrated  numerically  by  the  Runge-Kutta  method^  on  the  high  speed 
computers . 

The  initial  conditions  are  taken  at  the  shock  wave;  for  a  given  atmos¬ 
phere,  speed,  and  shock  inclination,  all  quantities  are  known  here.  The 
terminal  condition  is  =  0.  This  condition  makes  it  convenient  to  use  V 
51s  the  independent  variable  instead  of  cp. 

Figs.  II. 1  through  II. 9  contain  data  for  conical  flow,  for  which 
many  properties  of  interest  can  be  accurately  represented  as  functions  of 
the  parameter  sin  9^,  given  the  free  stream  conditions.  Pressure, 
temperature ,  density,  and  species  concentrations  on  the  body  surface  are 
plotted  against  M  sin  9^,  the  shock  wave  angle  is  also  presented  in  this 
form. 

5.  STAGNATION  VALUES  BEHIND  NORMAL  SHOCKS 

Calculations  are  made  of  the  flow  variables  when  the  fluid  is  brought  to 
rest  behind  a  normal  shock  wave,  as  for  instance,  at  the  intersection  of  the 
axial  streamline  with  the  surface  of  a  symmetric  blunt  body  (point  0  in 

Fig.  3*l*a).  The  independent  variable  in  eqs.  (2.6),  (2.9)?  and  (2.10)  is 
taken  to  be  the  velocity  q.  These  equations  plus  the  momentum  relation 

dp/dq  =  -  pq 

form  a  set  of  differential  equations  which  are  integrated  numerically  from 
q  =  q^  to  q  =  0,  the  terminal  values  of  the  variables  giving  the  stagnation 

conditions . 

Figs.  III.l  through  III. 9  present  the  thermodynamic  variables  and 
species  concentrations  for  stagnation  flow  behind  normal  shock  waves  as 
functions  of  M^  .  This  information  is  useful  in  studying  the  subsonic 
region  between  the  surface  of  a  two-dimensional  or  axisymmetric  blunt  body 
and  the  detached  shock  wave  ahead  of  it. 

6.  GRADIENTS  AT  SHOCK  ON  STAGNATION  STREAMLINE 

The  flow  variable  gradients  along  the  central  streamline  behind  a 
curved  shock  (point  P  in  Fig.  3.1. a)  can  be  calculated  if  the  curvature  of 
the  shock  is  known.  If  ct  is  arc  length  along  the  shock  wave,  it  is  seen 
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from  the  relation 


h  -  Kv  Ik  =  [cos(B  -  9)i  si  +  [sln(B  -  e)]  h 

(where  s  and  n  are  arc  lengths  along  streamlines  and  their  orthogonal 
trajectory,  respectively)  that 

W«0B  =g0°  •  (6.1) 

The  momentum  conservation  equation* 

pq2  Bfi/ds  =  -  dp/dn  (6.2) 

shows  that  (dp/dB)g  _  =  0.  Then,  by  Eq.  (2.6)  and  Eq.  (3.4)  differen¬ 

tiated  with  respect  to  8,  dp/dB  =  dT/dJ3  =  dc^/dB  =  dc^/dB  =  0  at  the  x-axis. 
By  the  differentiated  Eq.  (3.1)  and  Eq.  (6.l) 

(se/an)p  =  90o  =  [1  -  p/P##]Kw  .  (6.3) 

Expanding  [(sin  8)/y]  near  the  x-axis,  noting  that  dy  =  do  sin  8  , 

(sin  e)/y  =  [d0/dy)  ^  Q  y  +  •  •  •  ]/y  =  (de/do)  =  0  +  ‘ 

Jw  Jw  Jw 

Therefore 

(sin  9)/yw  5;Kw(de/d8)B  =  ^o  =  [1  -  (p/P<J!))  ]KW  *  (6.4) 


Substituting  Eqs.  (6.3)  and  (6.4)  into  the  flow  equation 

I  $£  _i_  *2  +  +  P  sin  9  _  o 

p  ^s  2  fts  an  y 

P<1 


(6.5) 


where  e  =  0  and  1  for  two-dimensional  and  axisymmetric  flow,  respectively, 
one  obtains 

dp/ds  =  q^  dp/ds  -  (l  +  e)K^pq2  [(p/p^  -  l)]  •  (6.6) 

On  substituting  dp/ds  from  Eq.  (6.6)  into  Eq.  (2.6)  the  gradients  of  the 
flow  variables  are  determined.  It  is  seen  that  the  gradients  are  all 
proportional  to  K  ,  and  that  for  a  given  M  the  axisymmetric  gradients  are 


*Eqs.  (6.2)  and  (6.5)  expressing  conservation  of  mass  and  momentum  are 
found,  e.g.,  in  Section  3  of  Kef.  5. 
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twice  those  for  two-dimensional  flow. 


Gradients  along  the  stagnation  streamline  at  the  shock  wave  are  shown 

in  Figs.  IV. 1  through  IV. 8  for  the  thermodynamic  variables  and  the  species 

concentrations.  Arc  length  is  given  in  terms  of  the  radius  of  curvature  of 

the  shock  wave,  R  ,  at  the  x-axis. 
w 

T.  COMPUTATIONAL  RESULTS 

Computational  results  are  presented  in  the  diagrams  which  follow. 

Graphs  of  desired  quantities  are  plotted  for  free  stream  temperatures  of 
273<l6°K  and  300°K,  and  free  stream  pressures  of  1.0,  .1,  .01,  .001,  and 
.0001  atmospheres;  the  choice  of  Mach  numbers,  cone  angles  and  shock  wave 
angles  makes  possible  a  temperature  range  coverage  of  35000°K  to  10,000°K. 

A  complete  survey  for  air  in  dissociation  equilibrium  is  not  feasible  be¬ 
cause  of  the  multitude  of  combinations  of  parameters.  It  is  felt,  never¬ 
theless,  that  the  following  set  of  diagrams  can  be  used  to  obtain  approxi¬ 
mate  information  adequate  for  planning  experiments,  predicting  and  checking 
experimental  results  over  a  wide  range  of  conditions  attainable  in  the 
laboratory. 

Comparisons  were  made  of  the  present  results  with  previously  published 
results  in  which  more  accurate  models  of  high  temperature  air  were  assumed. 

For  example,  conical  flow  parameters  (Fig.  II. 1  -  -  II. 4)  practically  coincide 
with  those  of  Romig  (Ref.  4).  A  comparison  of  shock  and  stagnation  pressure 
calculations  with  those  of  Feldman  (Ref.  l)  is  shown  in  Fig.  7.1;  the  largest 
discrepancy  of  all  the  parameters  is  found  in  the  stagnation  pressure. 
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FIG.  3.1.  CROSS-SECTION  DIAGRAM  OF  FLOW  FIELD 
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FIG.  1.1.  VARIATION  OF  PRESSURE  ACROSS  SHOCK  WAVES 
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FIG.  1.2.  VARIATION  OF  DENSITY  ACROSS  SHOCK  WAVES 
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FIG.  1.11.  VARIATION  OF  NO  CONCENTRATION  ACROSS  SHOCK  WAVES 
(Wco  =  28.86  GM/MOLE) 
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FIG.  IV. 6.  GRADIENT  OF  0  CONCENTRATION  ALONG  STREAMLINE  NORMAL  TO 
CURVED  SHOCK  WAVE  (U«  =  28.86  GM/MOLE) 
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FIG.  IV. 8.  GRADIENT  OF  NO  CONCENTRATION  ALONG  STREAMLINE  NORMAL  TO  CURVED 
SHOCK  WAVE  (Woo  =  28.86  GM/MOLE) 
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FIG.  7.1  COMPARISON  OF  PRESSURE 


CALCULATIONS  WITH  VALUES  IN  REF.  1. 


REFERENCES 


1.  Saul  Feldman,  "Hypersonic  Gas  Dynamic  Charts  for  Equilibrium  Air," 
Research  Rept.  40,  AVCO  Research  Lab  (1957). 

2.  Henry  E.  Hudgins,  Jr.,  "Supersonic  Flow  About  Right  Circular  Cones 

at  Zero  Yaw  in  Air  at  Thermodynamic  Equilibrium,"  T.  M.  1493,  Picatinny 
Arsenal,  Dover,  N.  J.  (1964). 

3.  Paul  V.  Marrone,  "Normal  Shock  Waves  in  Air:  Equilibrium  Composition 
and  Flow  Parameters  for  Velocities  from  26,000  to  50,000  Ft/Sec," 

CAL  Rept.  No.  AG-1729-A-2,  Cornell  Aeronaut.  Lab  (1962). 

4.  Mary  F.  Romig,  "Conical  Flow  Parameters  for  Air  in  Dissociation 
Equilibrium,"  Research  Rept.  7>  Convair  Scientific  Research  Lab.  (i960). 

5.  J.  Spurk,  N.  Gerber,  and  R.  Sedney,  "Characteristic  Calculation  of  Flow 
Fields  with  Chemical  Reactions,"  Aberdeen  Proving  Ground,  ERL  R-1276 
(March  1965).  Also  AIAA  J. ,  4,  50-57  (1966). 

6.  Paul  V.  Marrone,  "Inviscid,  Nonequilibrium  Flow  Behind  Bow  and  Normal 
Shock  Waves,  Part  I.  General  Analysis  and  Numerical  Examples,"  CAL 
Rept.  No.  QM-l626-A-12(l) ,  Cornell  Aeronaut.  Lab.  (1963). 

7.  K.  L.  Wray,  "Chemical  Kinetics  of  High  Temperature  Air,"  Progress  in 
Astronautics  and  Rocketry,  Vol.  7,  Hypersonic  Flow  Research, 

F.  R.  Riddell  (Editor),  Academic  Press,  New  York  (1962). 

8.  Ames  Research  Staff,  "Equations,  Tables,  and  Charts  for  Compressible 
Flow,"  NACA  R-1135  (1953). 

9.  A.  Ralston  and  H.  S.  Wilf,  Editors,  Mathematical  Methods  for  Digital 
Computers  (p.  110 ),  John  Wiley  and  Sons,  New  York  (i960). 


56 


APPENDIX 


COEFFICIENTS  OF  EQUATIONS  (2.6)  AND  (2.9) 


A11  = 

2/p 

A12 

A13  =  2/^C0  +  CN 

+  (2/W  )] 

CD 

ai4  = 

2/T 

c  r!o  +  Sli 

A21  " 

CoI_Kl  2K^  1 
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2pc0  i  PCN 

Kx  2  2K 

A23 

pco 

2K5 

A24  = 
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ICO 
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1  1 
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CL 
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CN  dK53 

2K52  dT  ' 

.  rcN  co  -] 

A32  = 

PCN 

2K5 

A3i 

Cn‘k2  2K5_ 

A33 

2PCN  +  1  +  PC0 

K2  2  2K5 

[CN  ^2  CQ 

A34  "  PCnL  2  dT  ov  2  dT 


A4l  ^C0  h02/Kl^  +  ^CN  +  ^CNC0  hNOy/K5'> 


TABLE  I 


REACTION  RATE  COEFFICIENTS  AND  EQUILIBRIUM  CONSTANTS 


Dimensions:  ( k  f )  j  --  m^/(kmol  sec)  ; 

Kj  for  3  body  reactions  --  kmol/m^ 
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TABLE  I  I 


SPECIFIC  ENTHALPIES  AND  SPECIFIC  HEATS 


V  7RT/2  •  <R602>/  CexplS^/TJ-l]  * 

hN2=  7rt/2  ♦  (R8„2)/  t=*P(SNg/T)-l]  • 

h0  *  6RT/2  *  *0 

h„  ■  5RT/2  *  hu 

N  N 

bN0=  7RT/2  *  (R8m0)/  [exp(6N0/T)-l]  *  hN£) 


<Cp)02=  7R/2  *  R(0Oo/T)2exp(eo!>/T)/[exp(9O2/T)-l]2 


(Cn)kl  =  7R/2  +  R(GN  /T)2exp(ew  /T)/[exp(0N^/T)-l]2 


P  N/ 


N< 


( C  p ) 0  =  5R/2 


(Cp)N  =  5R/2 

(Cp^N0=  7R/2  +  R*®N0^  exp(0NQ/T)/ [exp(0NO/T)- 1] 


No. { i ) 

Spec i es 

6j(°K) 

h  (Dyn  m/k  mole) 

1 

°2 

2256 

0 

2 

n2 

337U 

0 

3 

0 

-- 

g 

2.U67  65  x  10 

U 

N 

U.7I0  63  x  I08 

5 

NO 

2719 

0.898  655  x  I08 

5 

I  Dyn .  =  10  dynes 
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